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Abstract 

The random flux model (defined here as a model of lattice fermions hop- 
ping under the influence of maximally random link disorder) is analysed field 
theoretically. It is shown that the long range physics of the model is de- 
scribed by the supersymmetric version of a field theory that has been derived 
earlier in connection with lattice fermions subject to weak random hopping. 
More precisely, the field theory relevant for the behaviour of n-point correla- 
tion functions is of non-linear a model type, where the group GL(n|n) is the 
global invariant manifold. It is argued that the model universally describes 
the long range physics of random phase fermions and provides further evidence 
in favour of the existence of delocalised states in the middle of the band in 
two dimensions. The same formalism is applied to the study of non-Abelian 
generalisations of the random flux model, i.e. A'^-component fermions whose 
hopping is mediated by random U(A^) matrices. We discuss some physical 
applications of these models and argue that, for sufficiently large A^, the ex- 
istence of long range correlations in the band center (equivalent to metallic 
behaviour in the Abelian case) can be safely deduced from the RG analysis 
of the model. 

PACS numbers: 71.30.-hh,72.10.-d,05.45.+b 
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I. INTRODUCTION 



Quantum disordered systems are typically described in terms of Hamiltonians of the form 

H = Ho + V, (1) 

where models the underlying clean system and disorder is introduced via the randomly 
distributed Hermitian operator V. Sometimes, however, it is preferable to describe disor- 
dered systems differently, viz. by Hamiltonians where disorder is implemented in terms of 
unitary stochastic operators. A simple prototype example is the system of lattice electrons 
subject to a (perpendicular) strong random magnetic flux (see fig. |I|). An electron hop- 
ping from one site to the next acquires a field dependent phase which, if the field is strong 
enough^ becomes completely random. Under these conditions the system is described by 
the Hamiltonian, 

^ = -E4e^''^c„ (2) 

where is a summation over all nearest neighbour sites of a dimensional bipartite 

lattice, Cj is an electron field operator at site i, (pij a phase living on the link i ^ j, and 
maximum stochasticity implies that the phase (pij is drawn randomly from the full interval 
[0, 2tt] subject only to the Hermiticity requirement (pij = —(pji. Notice that the definition 
(0) implicitly assumes that all energies are measured in units of the modulus of the hopping 
matrix elements. 

It has been suggested that the model (Q) describes the physics of the fractional quantum 
Hall effect close to half filling (where the problem is expected to map onto a system of 
composite fermions subject to a strong random field with zero mean ||T|]). Furthermore, it 
is believed to be of some relevance in the gauge theory of high- Tc superconductivity 0. 
(However, it is important to recognise that the denotation 'random flux model' is, to some 
extent, misleading in the sense that the Hamiltonian (^ does not actually model a short 
range correlated random magnetic field but rather the lattice version of a random vector 
potential. The difference can be of importance since a short range correlated field would 
imply the presence of long range correlations in the vector potential Keeping this 

difference in mind we nevertheless continue to refer to (|^) as a random flux Hamiltonian.) 

The random phase Hamiltonian (^), describing the lattice dynamics of sca/ar fermions, 
can straightforwardly be generalized so as to accommodate fermions with an internal struc- 
ture: Let us assume that a) the lattice fermions carry an internal quantum number, 
Cj — > c",a = 1,...,A^ such as, e.g., spin and b) that the fundamental hopping mecha- 
nism leads to a full randomisation of the internal states. Under these circumstances, the 
Hamiltonian (^ generalises to (c.f. fig. ^ 

H = -J24U^,C,, (3) 



More precisely, maximum stochasticity requires fields stronger than one flux-quantum per lattice 
cell. 
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where Uij = {f/,"'^} G U(A^) are A^-dimensional unitary matrices sitting on the nearest neigh- 
bour hnks of the lattice, and q = {cf^}, a = 1, . . . , N, a = 1, . . . ,n are x ra-component 
field operators. The lower index a, which has been introduced for future reference, plays 
a passive role in the sense that it is summed over and does not couple to the f/'s. 'Full 
randomisation' upon each hopping process translates to the Uij being drawn democratically 
from the Haar measure of the group (subject to the Hermiticity requirement Uij = Uj^.) 
Following a convention introduced in Ref. we refer to the index a as a 'colour' index, 
whilst a will be denoted as a 'fiavour' index (for the motivation of this terminology, see 
below). The original RF-Hamiltonian fits into the general framework in the sense that it is 
the N = 1 version of @). For the sake of simplicity we will also refer to (^) as a 'random 
fiux' model. 

The generalized RF-model has applications in both condensed matter and high energy 
physics: 

Physics of manganese oxides: Suggestions have been put forward to explain the effect of 
giant magnetoresistance, shown by transition metal oxide compounds, in terms of effective 
models consisting of lattice electrons hopping in a matrix of localised 5 = 3/2 core spins [Q. 
The hopping process is essentially influenced by the Hund's rule coupling between the spin 
of the itinerant electrons and the localised spins. In the limit where the coupling strength 
becomes infinite, the hopping matrix element along any given bond i ^ j depends solely on 
the orientation of the spins at the sites i and j 0, i.e. the hopping is mediated by an SU(2) 
matrix which acts in the 5 = 1/2 space of the electron spins and depends on the orientation 
of the two adjacent core spins. Assuming that a) the spin orientations are random, b) that 
the effect of the randomness can be mimicked by drawing the link matrices from a random 
distribution, and c) that, with regard to universal properties, the difference between random 
SU(2) hopping and random U(2) hopping is of no concern, we arrive at the (A^ = 2,n = 1) 
version of (|). 

Lattice gauge theory: The simplest^ lattice action of a matter field with A^ 'colour' indices 
subject to a non-Abelian local gauge symmetry is given by 



(i,j) ^9 □ 



ltro(f/^))+h.c. 



(4) 



Here, a = 1, . . . ,n has the physical significance of a fermionic fiavour index, X^n is a sum 
over all plaquettes of the lattice, and tr^ ([/■*) = tr {UijUjkUuUu) represents the trace over 
all f/'s sitting on the links of the plaquette. Physically, the plaquette action represents the 
discrete version of a continuum field strength tensor (c.f., for example, |^). In the strong 
coupling limit, g ^ oo (i.e. the long distance-confinement phase of the theory) the gauge 
fields become free and (HI) is governed by a structure like (0). 



more realistic lattice representation of the QCD Lagrangian would involve Dirac 7-matrices 
acting in the Dirac-spinor space of the fermionic matter field. Although the embedding of this 
additional structure into the formalism discussed below is straightforward, we prefer to keep the 
discussion simple and not to include the Dirac-spinor structure. 
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What can be said, in general terms, about the physical behaviour of the scalar random 
flux Hamiltonian and its non-abelian generalizations?: 



At flrst sight it looks like the random flux Hamiltonian falls into the general category 
of Anderson Hamiltonians for disordered systems. More precisely, the model (|^) bears 
similarity with a lattice Anderson Hamiltonian with strong disorder in the hopping 
matrix elements. Being Hermitian but not real, the Hamiltonian (Q) is expected to 
fall into the general class of models with unitary symmetry. In particular, all states 
might be expected to be Anderson localised. As will be discussed below, this picture 
is correct as long as one stays away from zero energy: 

The middle of the tight binding band, e = 0. In the band center, the Hamiltonian 
possesses a discrete symmetry which, as will be explored in detail below, heavily 
influences the physical properties of the model. The nature of the symmetry can easily 
be understood by inspection of flg. |^: The essential point is that the Hamiltonian of 
the RF problem is of pure nearest neighbour type. (There are no matrix elements 
connecting any one site with itself.) To connect this statement with the existence of 
a symmetry, one may subdivide the lattice into two sublattices (the o's and the x's 
in fig. When represented in a o/x block decomposition the Hamiltonian becomes 
purely off- diagonal, 

H=(fr (5) 



H 



X,o 



Formally, the block 'off-diagonality' of the the Hamiltonian can be expressed as 

a^Has = -H^ [H, a-^U = 0, (6) 



where [ , ]+ is the anticommutator and as a Pauli matrix in o/x space. Symmetries 
of the type (§) are commonly referred to as chiral symmetries. 

The existence of a chiral symmetry has far reaching consequences for the physics of 
the model at the band-center. Before turning to a more comprehensive discussion of 
this point, let us summarise a few elementary facts that follow straightforwardly from 
the definition of the model. Firstly, we note that the Hamiltonian is always chiral, 
i.e. not just in the band center. One may then ask why, for finite energies, the chiral 
symmetry is ineffective. The reason is that the information about both thermodynamic 
and transport properties of the system is encoded in the single-particle Green function 
G^{e) = (e^ — H)~^ rather than in the Hamiltonian itself (here = e ± iO). Now, 



i.e. the inversion property under adjungation with is broken by finite e. Some 
physical implications of this symmetry breaking will be discussed in detail below. 

Secondly, we note that stochastic Hamiltonians with a chiral structure have previously 
been studied in the random matrix literature P|-p^. It has been shown that they 
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belong to symmetry classes that are fundamentally different from the standard Wigner 
Dyson classes. In particular, the band center RF model belongs to the universality class 
of the so-called 'chiral GUE', ChGUE [jgiHH^ . For finite e, the model undergoes 
a symmetry crossover from ChGUE to ordinary GUE symmetry. 

To the best of our knowledge the full extent of the consequences of the chirality of the 
RF model was first noticed in a seminal paper by Miller and Wang ||16|| . 



Beginning with Lee and Fisher ^7j, there has been intensive research into the An- 
derson localization behaviour of the model. So far the question whether or not all 
states are localised has not been conclusively settled: Despite extensive numerical and 
analytic studies, opinion is divided between those who conclude that all states are 



localised |18,19,4,20-22|, and those who believe that delocalised states persist at the 



center of the band [p^- ^ , p!6| , ^9| -0ll] . Below we will provide further evidence for the ex- 



istence of a metallic phase in the middle of the band. Needless to say, the phenomenon 
is induced by the unconventional symmetry of the model. By contrast, for finite e, the 
model can be shown to map onto the standard Anderson model of unitary symmetry 
which is known to be localised in dimensions two and less. For weak random magnetic 
fields, this connection has already been established microscopically in Ref. |Q. 

Besides the localization behaviour, several other physical properties of the RF system 
are affected by the chiral symmetry. In particular the presence of the symmetry gives 
rise to non-trivial quantum-interference corrections to average sm^/e-particle properties 
such as, e.g. the density of states. Some of these features will be discussed below. 

The subject of this paper is the construction of an effective field theory describing the 
low energy/long range properties of the generalized random flux model (^. Some empha- 
sis will be placed on exploring the consequences of the chiral symmetry of these systems. 
Specifically, it will be shown that the random flux model maps onto a field theory of a 



structure similar to that constructed earlier by Gade |^2[ in connection with weakly disor- 
dered sublattice models introduced by Oppermann and Wegner [^. Both the symmetries 
and the structure of Gade's field theory are known to be essentially influenced by chirality. 
In particular, renormalization group (RG) flow of the theory is drastically different from 
the one observed in conventional low dimensional disordered systems. Indeed, the fact that 
the RF models @ and the Wegner/ Oppermann sublattice systems of broken time reversal 
invariance, respectively, map onto the same field theory might have been anticipated from 
the outset, the reason being that the latter represents the maximally disordered counterpart 
of the former. On the field theoretical side, the main difference between the models lies in 
the value of the coupling constants. The consequences of this difference for the large scale 
analysis of the RF model will be discussed below, after the field theory has been microscop- 
ically derived. A second difference to Gade's analysis is that the present construction will 
be based on the supersymmetry formalism (as opposed to Gade's boson replicas). As will 
be discussed below, the supersymmetric formulation enables one to explicitly demonstrate 
that the ergodic limit of the model is described by chiral random matrix theory. 

The paper is organized as follows: We start out from a qualitative explanation for the 
importance of the chiral symmetry in section |I|. The central sections |T| and |V| contain the 
derivation of the long range field theory of the model. In section |V| we discuss what kind 
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of physical predictions follow from the field theoretical modelling. To clarify our discussion, 
much of the technical content of the analysis has been moved to appendices. 



II. THE CHIRAL SYMMETRY: QUALITATIVE CONSIDERATIONS 

Prior to embarking on the quantitative analysis of the RF model, it is instructive to 
develop some intuition regarding the consequences of the presence of the chiral symmetry. 
Perhaps the most apparent manifestation of the symmetry is the non-triviality of ensemble 
averaged single-particle properties such as the DoS: As shown in Ref. ||32| within a pertur- 
bative renormalization group scheme, the DoS of the two-dimensional (weakly disordered) 
sublattice model diverges upon approaching the middle of the band. The same type of be- 



haviour has been predicted for one-dimensional chiral systems 1^^. That these results can, 
in fact, not hold true all the way down to e = follows from exact random matrix theory 0. 
The random matrix analysis which, by general symmetry arguments, is expected to effec- 
tively describe the ergodic limit of chiral problems (i.e. the asymptotic vicinity of the band 
center) yields a vanishing DOS at e = 0. Since the breakdown of the DoS takes place on an 
energy scale set by the level spacing, it is generally difficult to detect within perturbative 
approaches. In fact it is known, for example from the conceptually related analysis of the 
level-spacing gap in the DoS of superconductor/normal metal systems without proximity 
effect |^5|, that a non-perturbative level spacing gap may signal its presence in terms of 
spurious divergences appearing in perturbative analyses. Thus, judging from perturbation 
theory alone, it may be difficult to tell what part of the band center DoS divergence is 'real' 
and what is merely an unphysical repercussion of a tiny non-perturbative effect. 

As will be discussed in more detail below, the complex behaviour of the DoS finds its 
origin in the fact that, in chiral models, mechanisms of large scale quantum interference — 
which are commonly represented by diffusion-type modes like 'diffusons' and 'Cooperons' 
— already operate on the level of single-particle properties. This should be contrasted with 
the situation in conventional disordered models where the average Green function picks up 
a largely structureless imaginary self energy but otherwise is not affected by the disorder. 

To understand why the DoS of the RF model behaves unusually in the band center, it 
is not even necessary to invoke the concept of diffusion modes. A simple path counting 



argument analogous to that employed by Brinkman and Rice [36| is sufficient: Consider the 



average DoS per unit volume, p, represented in terms of the single-particle Green function. 





' e- -H 

where the outer brackets represent the phase averaging 

1 r^'^ 

(id) ' 

Expanding in powers of H and inserting lattice resolutions of unity, J2j we obtain 



(7) 



(■■■>-n^f (8) 
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where L is the hnear extension of the lattice and a the lattice spacing. The expression 
Em=o '/'im.im+i has the meaning of the total phase acquired during the propagation along 
the path i ^ ii ^ 12 ^ ■ ■ ■ ^ i (c.f. fig. ^. The second equality in states that the 
contribution of all paths with a non-vanishing total phase averages to zero. The surviving 0- 
phase paths are self retracing — indicated by the symbol 's.r.'. In the third equality we have 
used the fact that the total number of self retracing paths of length n is 

(2rf)"/2. Evidently, 

the summation over n in the last line diverges for e < \/2d, indicating a breakdown of 
perturbation theory in the vicinity of the band center. Notice that the presence of 0- 
phase paths is intimately connected to the sublattice structure of the system. No analogous 
contribution to the mean DoS of conventional (on-site) disordered systems exists. 

A second and more formal way of demonstrating the non-triviality of single-particle 
properties is based on a few algebraic manipulations of the single-particle Green function. 
Introducing a system of Pauli matrices in the two-component 0/ x-partition, the Hamiltonian 
takes the form 



H 



Hi~ iH2 
Hi + iH2 




where U = exp(z^cr2) and both Hi and H2 are Hermitian. Similarly, the Green function can 
be represented as 



e+ - Hi +1H2 
-iH2 e+ + Hi 




The point now is that the central matrix can be interpreted as a retarded and an advanced 
stochastic Green function g~^{e) = {e'^ — Hi)^^ and g^{—t) = (— e+ — Hi)^^, coupled by the 
operator iH2. Since the the manipulations with the matrices U and do not affect the 
structure of the singularities of the Green function, we conclude that the single-particle G~^{e) 
has the effective complexity of two coupled stochastic Green functions. In particular, we can 
anticipate that the average G~^{e) will be affected by mechanisms of quantum interference 
between ^''''(e) and g~{—e). Notice that both large energy arguments e and on-site matrix 
elements of the Hamiltonian (corresponding to diagonal blocks in the 0/ x-partition) would 
lead to a destruction of this two-particle coherence. 

After these general considerations, we proceed to the quantitative study of the generalized 
RF-problem. Our general strategy parallels that employed by Zirnbauer [^] in a loosely 
related study of the mapping of the Chalker-Coddingtom network model of the Integer 
Quantum Hall Effect onto the effective field theory proposed by Pruisken. 
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III. MAPPING THE RANDOM FLUX MODEL TO A FIELD THEORY: 
COLOUR-FLAVOUR TRANSFORMATION 



As a first step towards the construction of a field theory of the RF problem, we need 
to introduce some degrees of freedom capable of describing the long range features of the 
model @. A mathematical tool for mapping model systems governed by unitary stochastic 
operators onto functional integrals containing a 'low energy' sector describing the large 
distance physics has recently been introduced and applied by Zirnbauer 0. For reasons 
that will become clear below, the transformation has been christened the 'colour- flavour 
transformation'. Our application of the colour- flavour transformation to the U(A^) RF- 
model closely parallels the analysis of Ref. of the spectral properties of single U(A^) 
matrices. For this reason the presentation of the initial stages of the construction of the 
theory will be comparatively concise but nevertheless self contained. (For a deeper discussion 
of the mathematics of the colour-flavour transformation we refer the reader to the original 
reference.) 

As a first step towards the transformed theory, we consider the Gaussian functional 
integral 

Z[J] = / / V{^,^) e <«^.^es> \ , (10) 

where 

(id) 

is the non-Abelian generalisation of the phase average (§), dUij denotes the Haar measure 
on U(A^) and A (B) is the sublattice of o-sites ( x -sites). The fields 

carry three types of indices: A lattice index i, a 'colour' index a (coupling to the internal 
group indices of the U^s), and a two component 'flavour' index a = {a, m). Here m = 1, . . . , n 
is a discrete counting variable whose significance will become clear shortly, whereas cr = 1, 2 
introduces a boson/fermion grading into the theory: Field components ifj^^^^^ = S°'^ 
are complex commuting whilst i'"„=2m = Xfm '^^^ elements of a Grassmann algebra 
('anticommuting variables'). Convergence of the integral ( p!OD requires that the components 
Sf^ be the complex conjugates of S"^, while the variables Xim ^'^d x^m independent of 
each other. In cases where no confusion can arise we will employ the more explicit notation, 
cr = B, F instead of o" = 1, 2. 

The elements of the diagonal matrix 

z = diag(2;i, . . . , 2;„) O ILsat, 

where (see below) G C, have the significance of energy arguments of single-particle Green 
functions. The presence of the matrix 
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s = diag(sgnlmzi, . . . ,sgnlmz„) ® I2N 

ensures convergence of the Gaussian integral. Finally, the 'source field' J is block diagonal 
in flavour space, 

J = bdiag (Ji, . . . , J„) (g) 

where Jm = {{jm)ij}, m = 1, . . . ,n are arbitrary matrices on the lattice, and is the Pauli 
matrix in the space of grading or a-indices. 

The signiflcance of the functional integral (p!0| ) lies in the fact that averaged products 
of single-particle Green functions can be obtained from it by differentiation with respect 
to matrix elements of the source fleld J: As with the application of supersymmetric fleld 
integrals to the analysis of Green functions of conventional disordered systems, we have the 
identity 

^"^^^}n = 2" ( ri G,^.^i^j) . (11) 

0{Jl)^^h■■■0{Jn)i„j„ \m=l / 

In the following it is important to notice the existence of a symmetry of the functional 
(|lO|) which is directly inherited from the chirality of the Hamiltonian: Consider a linear 
transformation of the integration variables, 

^j^i^jV^ J^B, (12) 

where T,T,V,V G Gl(n|n) are arbitrary invertible supermatrices. In the limit z 0, the 
(source free, J = 0) action remains invariant provided that the transformation matrices are 
related to each other by 

TsV = s, VsT = s. (13) 

In other words, the z = functional integral possesses a global Gl(n|n) x GL{n\n) sym- 
metry. Notice that for a general non-chiral Hamiltonian (i.e. a Hamiltonian including 
diagonal blocks in the representation (^), the symmetry of the functional integral collapses 
to GL(n|n). 

What is missing in the above argument is that the matrices T, T, V and V can in fact 
not be chosen independently: The relation ip^ = implies that Vbb = ^bb = 

Tgg. Consequently, in the BB-sector of the theory, the global symmetry collapses from 
GL(n) X GL(n) to just GL(?t,). Later on we will see that the symmetry manifold is even 
smaller than that. We will come back to discussing the consequences of such reduction 
mechanisms after the theory has been subjected to the colour-flavour transformation. 

Having identifled the basic symmetries of the functional integral, the next conceptual 
step in the construction of the fleld theory is the average over the bond matrices Uij. At this 
stage the 'colour-flavour transformation', mentioned above, is applied. It effectively trades 
the group integrals over the Uij for integrations over a set of auxiliary fleld variables Zij: 
Applying the identity Eq. (3) in Ref. 0, i.e. the transformation identity for a single U(iV) 
integration, to all links of the lattice we obtain 
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Z[J] = J^iZ, Z)P(7/',^/')e^^<»eA,,6s> 



strln{l-ZijZij) 



xe «eA V « / V J / • _ (^-^^^ 

Here 

Zij = {{Zij)aa'} G GL(?2|n), Zij = {{Zij)aa'} G GL{n\n) (15) 

are 2n x 2n-dimensional supermatrices living on the directed link from the A site i to the 
nearest neighbour 5-site j. (By 'directed' we mean that there is no Zji.) The symbol 
J2jeNi stands for a summation over all j's that are connected with i and 'str' is the standard 
sup ert race |]38[| . 

In (0), the integration JV{Z,Z) over all pairs {Zij,Zij) is subject to a number of 
constraints: The fermion-fermion blocks of the matrices Zij and Zij are anti-Hermitian 
adjoint to each other, 

^ii,FF = ~^ij-,FF- (16) 

In the boson-boson block, a similar but more subtle constraint needs to be imposed. In the 
original formulation of the colour-flavour transformation the BB-blocks of Z^j and Zij 
were chosen to be Hermitian adjoint, Zij^BB,orig. = -^ij-.BB.orig.- will become clear later on, 
this convention is not suitable for the analysis of the RF-problem. We rather need to chose 
a different integration manifold, where 

Zij, BE = Zjj^^^ =^ Zjj^^^. (17) 

In order to define the meaning of the symbol we first note that Zbb G GL{n). Further- 
more, due to GL(n) = Gl(n)/U(n) x U(n), any matrix Z^b can be decomposed as 

Zbb = e^^'e^, (18) 

where both H and H' are Hermitian (implying that the factor G U(n)\GL(?7,) ~ 
GL{n)/\J{n) whilst e*^ G U{n)). The |-operation is then defined by 

Zt = (e^^^'e^)t ^ e-^e^^', (19) 

i.e. it changes the signs of the exponents in a manner that is opposite from the standard 
Hermitian adjunction 'f. As detailed in Appendix |^, the new integration domain can 
be reached from the original one by analytic continuation. Hence, switching between the 
two conventions does not affect the result of the integral. Finally, the Grassmann valued 
components, Zij -^F, ^ij,FB, ^ij,BF, ^ii.FB are independent n x n-matrices. 



Conceptually, the transformation has a status similar to the Hubbard- Stratonovich 
transformation employed in the standard construction of field theories for systems with 
Hermitian disorder ||39|: The motivation for introducing the Z- fields is that they couple 



to the fields '0 in a way different from the coupling ijjUijj between the ^p^s and the unitary 
disorder (see fig. §). Qualitatively, the Z's connect pairs of fields that, by construction. 
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represent segments of 'paired paths'. As discussed in section |T|, segments of this structure 
are stabihsed by mechanisms of quantum interference, implying that the Z's connect to a 
sector of the theory which contains information on long range behaviour. In addition, the 
Z's couple non-trivially to the 'flavour' space thereby potentially allowing for the description 
of quantum interference between different flavour species different Green functions). The 
subsequent analysis of the field theory will essentially amount to extracting the low energy 
interference stabilised sector contained in the Z-functional. 

The first step towards an effective long range model is to integrate out the ^/'-fields. At 
this stage one runs into an unexpected problem, viz. the Gaussian integration over the 
bosonic component of the ip^s does not converge for general complex energy arguments Zm- 
In other words, it looks like the orders of integrations over Z and ip cannot be exchanged 
which would imply that the colour-flavour transformation is useless. 

The problem can be surmounted by analytically continuing the arguments from their 
natural definition range = =t i^, ^ positive infinitesimal, a real 'energy' argument, 
deeply into the complex plane: em^iS e±2id. As shown in Appendix |B|, this manipulation 
cures the problems associated with the ^/'-integration. Moreover, we will see shortly that the 
remaining integration over Z is not plagued by convergence problems (as long as the Zm^s 
stay away from the real axis.) As a consequence the shift into the complex plane can be 
redone immediately after the ip^s have been integrated out. 

Integrating the partition function ([T^) over the field (the manipulation of the energy 
arguments mentioned above being understood), we obtain the functional integral 

Z[J] = jv{Z, Z)e-^[^'^lX[J], 

S[Z, Z] = -N J2 str In (1 - %%) + 

+Ar^strln lz+ zA +Ar^strln [z + ^iA , (20) 

where X[J] represents certain pre-exponential contributions resulting from an expansion in 
powers of the source field J to nth order. Since, in this study, the focus will be on an analysis 
of the general structure of the theory rather than on the computation of specific correlation 
functions, we will henceforth ignore this contribution. (To actually compute correlation 
functions one may proceed in a standard manner, i.e. keep track of -^[^], differentiate with 
respect to J at any convenient stage, and eventually carry out the Z-integral weighted with 
the effective action derived below.) 



IV. EFFECTIVE FIELD THEORY 



So far all manipulations have been exact. In order to make further progress towards an 
effective long range theory, we next subject the functional ( pUD to a standard saddle point 
analysis, that is we seek for solutions of the equations 



5S[Z, Z] 
5Zij 



0, 



Z=Zo,Z=Zo 
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5S[Z, Z] 



= 0. 

Z=Zo,Z=Zo 



Anticipating that spatial fluctuations of the fields Z and Z will be energetically penalised, 
we make an ansatz where the saddle point fields Zq and Zq connected to the A and the 
B sublattice, respectively, are spatially constant. Furthermore, the fact that the action is 
diagonal in the composite index a = (m, cr), permits one to seek for solutions of the form 

Substituting this Ansatz into the functional integral and varying the action we obtain 
the equations 

~Z, = Z, = -^z±x[-l + ^z'^ ' , (21) 

where x = {2d — 1)~^/^. Close to the middle of the band, 1, the solutions of (^iD 

compatible with the analyticity requirements imposed by the logarithms in the action are 
given by 

Zo = Zo^ ixs. (22) 

As in the construction of the field theory of the weakly disordered sublattice model P^ , 
the diagonal matrices (p^) are not the only solutions of the saddle point equations. The 
isotropy of the z = action in the space of a-indices rather implies that the set of solutions 
is highly degenerate. In order to identify the structure of the entire saddle point manifold we 
proceed in a standard manner and investigate how a general chiral symmetry transformation 
(p^) affects the diagonal saddle point solutions. From the structure of the action ([T^) it 
follows that under (|^) 

Z, = ixs-^'h'^'^ ^ ixs-^'^fTs~^/\ 

Z, = ixs~^'h-^'^ ^ ixs-^'^VVs-^'^ ixs^'\fT)-h^'\ (23) 

Since the action of the theory is invariant under global chiral symmetry transformations, all 
pairs (Zq, Zq) of the structure ( p3D are solutions of the saddle point equations, too. Notice 
that in ( p3D the Gl(r;,|n)-matrices T and T appear in a product combination. Taking into 
account that G'L{n\n) ■ GL(n|n) = G\{n\n), we conclude that the full extent of the saddle 
point manifold is given by GL(n|r7,).0 



^This is an example of what in QCD is known as the phenomenon of 'chiral symmetry breaking': 
Neglecting quark masses, the microscopic Yang-Mills action (corresponding to our -(/'-action in the 
small z-limit), is invariant under G x G where G is a symmetry group whose detailed structure 
depends on the flavour content of the theory etc. (In our formalism, its role is played by GL(n|n).) 
Now, in the effective low energy QCD Lagrangian, the microscopic chiral symmetry is spontaneously 
broken down to G x G ^ G, very much like the symmetry group of our low energy Z-functional is 
GL(n|n) rather than GL(n|n) x GL(n|n). 



12 



Eq. ( ]23|) defines the maximum manifold of solutions of the saddle point equations. In 
order to complete the determination of the global structure of the theory, we next need 
to explore how the manifold (^) intersects with the domain of the field integration (as 
specified by ([T5|) , (|16D and (|17D). To this end, we represent a general pair of matrices Z and 
Z according to 

{Z,Z) = ix{PT,T-^P), (24) 

where P,T & GL{n\n) and we have omitted the site index (ij) for notational transparency. 

In order to give the reparameterisation (^4]) some meaning, we a) have to be specific 
about the internal structure of the factor matrices T and P, and b) tell how the integration 
measure T^{Z,Z) is affected by the transformation. As for a), the restriction (|T^) leads to 
the relation 

Pff e GL(n)/U(n), Tff G U(n). (25) 

Similarly, the decompositions (0) and (|18|) imply that 

Pbb e U(n), Tbb e GL(n)/U(n). (26) 

The off-diagonal blocks Pbf^-Pfb^^bf and TpB are independent Grassmann variables. Fi- 
nally, counting the number of independent parameters, we see that the (P, T) manifold 
covers the entire {Z, Z) domain. 

As for the question of the measure, first note that Eqs. (p3D and (p^ imply that the 
saddle point configurations of the theory are specified by P = U. Consequently, excitations 
in P are massive, and can be treated in a quadratic approximation. Due to the fact that 
close to P = U the integration measure is approximately flat, we are spared addressing point 
b) explicitly. More precisely, the measure of the P-integration is fiat Euclidean, whilst the 
T-integration is over the invariant measures on the manifolds ( P^D and (^), respectively. 

Eqs. (15), (p5D and (p6|) contain the answer to the question of the general structure of 
the target manifold of the field theory raised above: The maximal manifold of solutions of 
the saddle point equation is given by GL(n|?T,). In the BB-sector, the intersection of this 
manifold with the actual domain of the field integration is given by GL{n)/U{n) C GL{n), 
and in the FF-sector, by U(n) C GL(n). The Grassmann sectors are not restricted. We 



note that these global structures have been worked out earlier , albeit in a more general 
manner based on Cartan's complete classification of symmetric spaces. Apart from the fact 
that it is not obvious how that more abstract line of reasoning can operationally be adapted 
to the RF-model, the present treatment has the advantage of encompassing the massive 
sector of the theory (the P's). As we will see later on, the RF-problem differs from standard 
(non-chiral) field theories of disordered systems, in that the massive P's have a significant 
impact on the final structure of the functional integral. 



Finally, we note that the original boson-replica |^2[ and fermion-replica based ap- 
proaches chiral disordered model systems led to field theories over GL(P?T,)/U(Pn) and 
U(Pn), respectively. Here P — denotes the number of replicas. This field manifold can 
be interpreted as the replica analogue of the BB-sector of our supersymmetric construction. 

We next employ the saddle point decomposition (E^) and (EB) to reduce the exact func- 



tional (PPD to a simpler effective functional describing the long range behaviour of the model. 
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Technically, the procedure will be based on a continuum approximation in combination with 
a gradient expansion around the spatially constant saddle point manifold. It is important to 
realize that this evaluation scheme is, in fact, problematic: As compared to analogous treat- 
ments of weakly disordered models, the gradient expansion is not stabilised by the presence 
of the small parameter 'disorder strength/Fermi energy'. Apart from the external parame- 
ters Zm, all energy scales are rather of the same order: 0{1). With the non-Abelian versions 
of the RF- model, N > 1, the situation is less problematic, the reason being that the mean 
field analysis is stabilised by the parameter However, for the = 1 random magnetic 
flux model, the expansion around the mean field configurations is not controlled by a small 
parameter. In this situation all we can do is provide some independent physical evidence for 
the belief that the field theory constructed below remains valid all the way down to A^ = 1.0 
Some arguments to that effect will be given towards the end of the paper. 

Keeping these words of caution in mind, we now turn to the discussion of the continuum 
expansion of the theory. On general grounds, it may be anticipated that after integration 
over massive modes, the effective action takes the form 

S[P, T] Seff [T] = S,[T] + Sa[T] + S^[r], (27) 

where Sz[T] is the contribution due to finite energies z, Sq[T] the action associated with 
fluctuations of the Goldstone modes (the T's), and ^^[T] is a residual action induced by 
the interaction between massive and Goldstone modes, respectively. We next assume that 
the field configurations relevant to the long range physics are smooth (differentiable). The 
different contributions to the action can then be computed by a gradient expansion fol- 
lowed by a continuum limit. Albeit conceptually straightforward, the explicit formulation 
of this program happens to be somewhat cumbersome. For this reason, the derivation of 
the different contributions to the action has been deferred to Appendix 0. As a result we 
obtain 

5efr[T] = --J i-^stiidT-'dT) + i ^ str (i(T + T-^)) + 

+ ^ str (T-^dT) str (T-^dT)^ + S^[T], (28) 

where a represents the lattice spacing, = N/id, = C/8d, and C denotes a numerical 
constant defined in Eq. ( |C16| ). Further, the discrete lattice summations J2i have been 
converted to integrals over smooth field configurations T. Eq. (pSD represents the main 
result of this paper. Sh[T] represents a boundary action whose structure depends on the 
geometry of the lattice. For a rf- dimensional hypercubic lattice with Ni sites in direction-i, 



^We remark that this extrapolation has the same status as the applicability of the Pruisken action 
to the Integer Quantum Hall transition. Strictly, the field theory can be formally justified only in 
the limit of high Landau level index, or large conductance. Yet, studies |37| of the phenomenological 
Chalker-Coddington model, corresponding to a single Landau band, indicate that the same action 
controls the fixed point of the Hall transition when the conductance is 0{e^ /h). 
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^bm = ^ E n(-)^''^'"'^"strln(T(.iL„...,.,L,)). (29) 

si,...,Stj=0,l i=l 

Notice that the boundary action depends on the 'parity' of the site-numbers, i.e. on whether 
or not the number of sites in the different extensions of the lattice is even on odd. Both the 
imphcations of this fact and the topological content of the boundary action will be discussed 
below. 



V. ANALYSIS OF THE FIELD THEORY 

What kind of information can be extracted from the effective action (PB])? A compara- 
tively complete answer to that question may be given in the 'zero dimensional' case, i.e. the 
case of systems of linear extension L in a regime where all energy parameters are such that 
Zm <^ a'^/L'^. (Recall that all energies are measured in units of the modulus of the hopping 
matrix elements.) The denotation 'zero dimensional' is motivated by the fact that, as in the 
case of conventional field theories of disordered systems, the boundary conditions to be im- 
posed on dT imply that the kinetic energy associated with spatial fluctuations is quantised: 
For constant field configurations it vanishes, while the first non-constant configuration has 
a fluctuation energy of 0{No? / L"^), the 'Thouless energy' of the system]^. (Notice that in 
d = 2, the Thouless energy is small of 0{A), where A is the mean level spacing. Thus, 
unlike in weakly disordered 2d systems, there is no extended region separating the onset of 
the zero mode regime z Ec = o? / L"^ from the fundamentally quantum regime z ~ A.) For 
energies z^ <^ / LF' spatially fluctuating configurations are effectively frozen out and the 
action simplifies to 

So[T] = -^|postr {z{T + T'^)) + S^[Tl (30) 
where T G GL(n|n) is a single matrix (the zero mode), and 

- 2. J2 • (31) 

To illustrate how the zero mode action can be applied to the analysis of physical observ- 
ables, we consider the DoS (|^) in the zero dimensional limit. Applying the general recipe 
for the computation of Green functions outlined above to the source field J = a^^, the DoS 
is represented by 



Pie) 



= jRe y'rfre-^o[^]str((r + r-i)a3^0. (32) 



^ That c? jl? plays the role of the Thouless energy Ec ~ vpljL'^ = VpT/L"^ {vp: Fermi velocity, I: 
elastic mean free path, r: elastic scattering time) of weakly disordered systems can be seen from 
the following rough estimate: Remembering that all energies are measured in units of the hopping 
matrix element, vp ^ a in the middle of the band. Further, for strong disorder, r ~ 1 is of the 
order of the band width. This gives ~ a? /L"^. 
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The integral depends sensitively on the presence of the boundary term 5'[T]b which in turn 
depends on the parity of the lattice. If the number of sites is odd (i.e. the number of lattice 
sites along each and every direction is odd), the boundary term takes the form 

5b[T] = strlnT, (33) 

otherwise the term vanishes. In both cases the explicit computation of the zero-mode integral 
(^) for the particular case of n = 1 is detailed in Appendix 0. Explicit calculation shows 
that, for a lattice with an even number of sites, the DoS takes the universal form 

ttVo 
2 

while if the number of sites is odd 



P(e) = ^(jg(7r6po) + J?(7r6po)), (34) 



P(e) = ^^-y^ (jAf(7repo) - JAr-i(7repo) J7V+i(vrepo)) , (35) 

where Jn denotes the Bessel function of order A^. 

In fact, both of these results have been anticipated in the literature. Firstly, Eqs. ( |3^ 
and (B^) have been derived earlier within the framework of a supersymmetric non-linear 



cr-model for the chiral random matrix ensemble ChGUE |T5|,E2|, and within pure random 



matrix theory 0-12|. Indeed, both expressions have been discussed in their connection with 



chiral symmetry breaking of the QCD Dirac operator |9|- [ll| , ^ . Moreover, in the present 
context, the validity of Eq. (|3^) has been investigated numerically . 



In fact, the coincidence between the low energy RF- field theory and the non-linear cr- 
model for chiral random matrices is not restricted to the mean DoS but rather extends to 
arbitrary n-point correlation functions. Equivalently, it can be said that in the low energy 
regime the RF-model is described by chiral random matrix theory, a fact that could have 
been anticipated from symmetry arguments alone. 

Eq. ( P^D also contains the answer to the question of what happens to the DoS upon 
leaving the domain of perturbation theory (c.f. Eq. (^) and below): Ultimately (i.e. on a 
scale comparable with the mean level spacing) it drops down to zero. Notice, however, that 
the zero dimensional analysis cannot exclude the existence of divergent behaviour prevailing 
down to energy scales ~ Ec (see below). 

At the same time, Eq. p^ ) suggests that on energy scales e > pg ^, the DoS approaches 
the constant value Pq. This prediction must be taken with some care: In c? = 2 it is 
meaningless because the domain of applicability of the zero mode approximation is set by 
Po ^ itself and the present analysis has nothing to say about energies larger than that. For 
d > 3, Pq should have the significance of a 'mean' DoS of the random flux model. Although, 
indeed, Eq. (|3T|) agrees parametrically with approximate expressions for the band center 
DoS of large dimensional hypercubic lattices |0, it is not clear whether po has any true 
physical significance for the physics of the RF models. The reason is that (|3T| ) represents a 
bare parameter which can be, in principle, subject to substantial renormalization. As will 
be discussed below, the extent to which the parameter changes under the RG flow is, at 
present, not clear. 

Leaving the zero-mode regime and returning to the more complicated, truly extended 
case, we first note that the field theory described by (pHf ) is not entirely original but rather 
has a number of closely related precursors: 
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Gade's analysis of the weakly disordered sublattice model led to a boson replica version 
of the present model, i.e. a theory over fields T G GL{nR)/\J{nR), where i? — is 
the number of replicas. The action for these fields coincided with (pS]), save for the 
important difference that, due to the weakness of the disorder, the coupling constant 
b^^ was parametrically larger than one. 

In high energy physics the action ( ]28|) has been suggested on phenomenological 
and symmetry grounds as relevant for the effective description of partially quenched 
QCD [0. In that context, the base manifold is 4 + 1-dimensional whilst the fields 
T G f/(nj + 1|1), where Uf is the number of quark flavours. 

A similar fleld theory (albeit without a (tr(T^^9T))^ operator) was derived in connec- 
tion with a model of Dirac fermions subject to a random non-Abelian vector poten- 



tial 1^ . That the random Dirac model is closely related to the RF-problem considered 
here follows from the fact that, when put on a lattice, Dirac fermions are equivalent 
to free fermions subject to an average vr-flux per lattice cell. In other words, the Dirac 
model corresponds to a weak random flux model centered around an average vr-flux 
background. 

As for the fleld theoretical description, apart from being deflned on a different fleld 
manifold (SU'=(Ar)/SU(Ar), where SU'(Ar) is the complex extension of SU(Ar)g), the 
model differs from ( p8l) in that it contains a Wess-Zumino-Novikov-Witten (WZNW) 
operator (however, see Ref. |^). While the symmetries of the flelds G GL{n\n) would. 



in principle, admit a supersymmetric version of a WZWN operator, we believe that 
the presence of a WZNW action is induced by the Dirac nature of the clean spectrum 
implying that there is no room for such a contribution in our model. (Cf. the absence 
of a Pruisken term in the action of a weak disordered electron gas subject to a random, 
zero average magnetic fleld.) To check that conjecture explicitly, we have computed 
the action of an SU(2) instanton type conflguration deflned in the FF-sector of the 
theory. No imaginary contribution to the action was obtained, a signature for the 
absence of a WZWN contribution. For the application of a similar recipe to determine 



the coupling constant of topological contributions to the action, see Ref. |37 



Finally, a fermion-replica analogue of the action (^) was obtained in connection with a 
lattice model of a disordered ci-wave superconductor pT| , ^ , ^ . The connection to the 



RF-model follows from the fact that the clean limit of the ci-wave model is governed by 
the presence of Dirac type excitation spectra and can therefore be effectively described 
in terms of a model qualitatively analogous to that of 7r-flux lattice fermions. 

Various conclusions regarding the physical behaviour of the RF-model, most notably 



about its localization behaviour, can be directly inferred from Ref. ||3^. There it was shown 
that the conductance of the weakly disordered 2d model at the band center (which is essen- 
tially determined by the coupling constant b~^) did not change under one- loop perturbative 



^Notice that for T E SU^(iV)/SU(iV), det (T) = 1 ^ tr {T-^dT) = which explains the absence 
of the (tr(r~^5T))^ operator in this field theory. 
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renormalization. This observation suggests that a non-locahsed state might exist in the mid- 
dle of the band. Since the stabihty of the perturbative RG merely relies on the smallness of 
the parameters 6, <C 1, its results can be straightforwardly carried over to the N ^ 1 
non-Abelian RF model: The one-loop renormalization indicates that for iV 3> 1 the strongly 
disordered RF model exhibits metallic behaviour at the band center. 

It is natural to ask as to what extent this prediction survives the limit ^ 1. Actually, 
one may adopt an even more hesitant attitude and ask whether the basic structure of the 
action (^) is stable as A^ — *• 1 (c.f. comments made in Appendix y in connection with the 
gradient expansion). We believe that the answer to the second question is positive whilst 
the answer to the first one remains, to some extent, a matter of faith. 

Addressing the second issue first, we note that, by a standard argument of operator 
relevancy, the infrared continuum action should be dominated by operators with the least 
number of derivatives. Of these, there are only two invariant combinations that are compat- 
ible with the basic symmetries of the model, namely the two operators appearing in (PH]). 
(E.g. the presence of an operator like ~ sti {diT~^djT),i ^ j is ruled out by virtue of 
reflection symmetry, etc.). Thus, the basic structure of the action is dictated by symmetry 
alone. As for the coupling constants multiplying the operators of the action, the absence of 
dimensionful parameters in the RF problem fixes their value up to numerical factors about 
which we need not be concerned. 

Now, turning to the more serious question of the validity the RG results in the Abelian 
case, from the results presented here, there appears to be only one useful test of consistency: 
The existence of a metallic state for A^ = 1 with strong disorder would at least be compatible 
with the behaviour in the limiting cases: A^ = 1 and weak disorder, and A^ 3> 1 with arbitrary 
disorder. Moreover, this conclusion does also seem to be consistent with the most recent 
numerical investigations PT|. More concrete analytical evidence in favour of band center 



delocalisation does not seem available. 

Besides the conductance, the DoS of the extended model also behaves in an unusual 
manner: The analysis of Ref. ||3^ predicts divergent behaviour upon approaching the band 



center, where the detailed functional form of the divergence depends on the dimensionality 
of the system. In particular, in two dimensions the RG analysis predicts 



g-/t^-ln(|e|) 

P(e) ^ • 

The scale below which the |e|^^ divergence begins to dominate over the exponential factor is 
determined by the constant k. Unfortunately the dependence of that constant on the relevant 
parameters of the model can not be reliably extracted from the perturbative RG analysis, 
i.e. the only statement that can safely be made is that a divergence forms somewhere upon 
approaching the band center. Ultimately, the divergent behaviour is cut-off by the onset of 
the non-perturbative ergodic regime discussed above. A quantitative understanding of both 
the formation of the divergence and its truncation in the close vicinity of the band center 
would necessitate a combination of an RG analysis with a non-perturbative treatment of the 
zero mode contribution, a program that is well beyond the scope of the present paper. 

A third piece of information that can be extracted from Ref. is how finite energy 
arguments = Cm ± ^0 affect the behaviour of the model: In the spirit of the RG, finite 
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energies have the significance of relevant perturbations. Thus, the band center represents 
a critical manifold whilst, for finite z^, the model flows into a different universality class. 
More precisely, the band center field theory — a model of universality class AIII in the 
classification of Ref. or, more commonly, a model of ChGUEl-symmetry — develops 



into an ordinary GUE type model. This is accompanied by a reduction of the symmetry 
manifold from Gl{n\n) to the standard GUE manifold GL(n|n)/Gl(n_|_|n+) x GL(n_|n_), 
where + = n and n± is the number of energy arguments with positive/negative 
imaginary part. (Note that this conclusion is in accord with the analysis in Ref. Q where 
a continuum random field model far away from the band center was shown to be described 
by a supersymmetric unitary model.) 

In that context, one may ask how the degrees of freedom of the model of reduced, 
unitary symmetry, commonly denoted by Q |^9[, emerge from the fields T of the present 



model upon crossing over to the unitary regime. In fact the full RF model could have been 
formulated in terms of Q-matrices (albeit Q's of a symmetry different from the standard 
unitary symmetry) from the outset: I.e. by defining Q = TAT^^, where A = ^ s, 
T = exp{iWa2) and the operate in sublattice space. The matrices T are then related to 
our T's by the identity T = exp{iW). The price to pay for this more common Q-formulation 
of the theory would be a certain amount of redundancy: A description in terms of structures 
Q ~ TAT~^ is the natural one for problems with coset space symmetry (e.g. the unitary 
coset space GL{n\n) /GL{n^\nj^) x GL{n_\n^)). However, the present problem has group 
type symmetry (Gl(n|n)) implying that a modelling directly in terms of group degrees of 
freedom is more natural. 



VI. SUMMARY AND DISCUSSION 

Summarising, we have mapped both the Abelian {N = 1) and non-Abelian {N > 1) 
random flux model onto an effective continuum field theory. The model differs in only two 
respects from the field theory obtained earlier for band center lattice fermions subject to 
weak off-diagonal disorder: It is a) a based on the formalism of supersymmetry (as opposed 
to boson replicas), and b) governed by a different coupling constant ~ x C(l) (as 
opposed to the large value obtained in the weakly disordered case). 

The supersymmetric formulation of the theory allows one to rigorously establish con- 
tact with the formalism of chiral random matrix theory |^14], and to demonstrate that 



the ergodic sector of the model is described by a random matrix ensemble of appropriate 
symmetry. Apart from this point, however, the fact that the present model is formulated 



in terms of supersymmetry, whilst the earlier analysis utilised a replica formulation is 
arguably incidental. What is, perhaps, more relevant is that both a weakly and a maximally 
disordered model with sublattice symmetry are described in terms of the same effective 
field theory. The only difference lies in the value of the coupling constant b~^. That this 
constant takes a large/small value in the case of weak/strong disorder follows readily from 
its significance as the bare conductance of the model. Given that the two limiting cases 
(very weak/very strong disorder) lead to the same effective action (^8]) , it is tempting to 
conjecture that the present field theory universally describes the long range behaviour of 
disordered sublattice models, where the amount of disordering is reflected in the value of 
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the coupling constant b ^. 

Unfortunately, the small value of the coupling constant in the strongly disordered case 



prevents us from directly transferring the results of the perturbative RG analysis of Ref ||32 
to the RF case. E.g. we cannot exclude the existence of some fixed point separating the 
delocalised regime for weak disorder from a potentially localised regime of the fully random 
flux model. Thus, the present analysis is not able to reliably answer the long-standing 
question of the behaviour of the strongly disordered Abelian RF-model in the band center. 
What can be said with some confidence is that the non- Abelian versions of the RF-model do 
exhibit long range correlations delocalisation), this being a consequence of the fact that 
b^^ oc implying that the RG is stable. (A more detailed account of the application of the 
present formalism to the physics of systems describable in terms of non-Abelian RF-models 
(or relatives thereof) will be presented elsewhere.) 

Finally, we mention that the present analysis of the RF-model can be generalized so as 
to encompass other lattice problems with unitarily implemented disorder: Firstly, models 
away from the 'strong coupling limit' (models with a finite plaquette action like in (^) can 
be dealt with at the expense of introducing extra 'ghost components' in the ip-field |[50|| . 
Secondly, models where stochasticity is implemented in terms of subgroups of the U(A^), e.g. 
as in the QCD-Hamiltonian (SU(A^) C U(A^)) or the random Ising Hamiltonian {Z2 C U(l)), 
can be analysed in terms of suitably adapted versions of the colour-flavour transformation 
| 5l| . Some of these generalisations are currently under investigation. 

We would like to thank A. V. Andreev, S. Hikami, M. Janssen, V. Kravtsov, P. B. 
Littlewood, E. Miiller-Hartmann, A.Tsvelik, J.J.M. Verbaaschot, and M. R. Zirnbauer for 
instructive discussions. 



APPENDIX A: CHOICE OF THE BOSONIC INTEGRATION CONTOUR 

In the original formulation ^ of the colour flavour transformation, the integration over 
the bosonic block of the fields Z, 

f dZ{. . .) 

extended over a manifold Mb defined by 

Mb = {{Z,Z)\Z = Z\ \EV{ZZ)\<1}. (A1) 

Here |EV {ZZ)\ < 1 is a shorthand notation for the 'modulus of all eigenvalues of the matrix 
ZZ smaller than one'. (As we will be concerned solely with boson variables throughout, the 
subscript 'BB' in Z has been omitted.) 

The Hermitian connection ( |A1[ ) between Z and Z ultimately finds its origin in the inter- 
pretation of these fields as parameter sets of bosonic spin coherent states. (A comprehensive 
discussion of the role of Hermiticity in connection with boson spin coherent states can be 
found in Perelomov's book |^2|-) '^^^ enforced by criteria of mathematical consis- 

tency. In order to encompass the low-lying modes of the RF-problem it is in fact necessary 
to abandon this connection and to consider pairs {Z,Z = Z^ ^ Z^\ where the matrix Z'^ 
has been defined in Eq. (^). More precisely, we need to analytically continue to a different 
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integration manifold jVb on which the connection between Z and Z = is still one-to-one 
but different from Hermit ian adjungation. 

Apart from the Hermitian connection between Z and Z, Eq. ( pO] ) contains a second 
condition, viz. |EV(ZZ)| < 1. As opposed to Hermiticity, this criterion is indeed needed 
for mathematical reasons. The most straightforward way to see this is to observe that the 
BB-block of the logarithmic term appearing in the colour-flavour transformed action ([I^ ) is 



given by tr In (1 — ZZ^ The condition |EV < 1 ensures that the cut of the logarithm 

is not crossed upon integration over Z. In fact, the condition can be relaxed to 

EV(ZZ) ^]l,cx)]. 

To understand the connection between the integrations over A^b and A/b, respectively, we 
reconsider the decomposition 



GL {n)3Z = WV, W = ex.p{H) G GL (n)/U (n), V = exp{iH') e U 



where H and H' are n-dimensional Hermitian generators. The BB-restricted integration 
over a function f{Z, Z) can then be written as 

/ f{Z,Z)= f dW f dVf{WV,V-'W), 

GL(n)/U(n) U(n) 
EV(W)<1 

where we have used the fact that Z = Z'^ = iWV)'^ = V~^W, and that the eigenvalues of 
ZZ = depend only on the Hermitian factor matrix W. The symbols / dW and / dV 
denote an integration over the invariant measures on the coset space GL (n) /U (n) and the 
group U(n), respectively. The basic idea of the analytic continuation is now to 'exchange' 
the role of the compact U (n) and the non-compact GL (n)/U (n) submanifolds, respectively. 
To this end, notice that both W and V can be unitarily diagonalised, 

W = exp{H) = Uexp{D)U^ ^ dW = dUdDJ{D), 

V = exp{iH') = U' exp{iD')U'^ ^ dV = dU'dD' J'{D'). 

Here U,U' G U (n) are the diagonalizing matrices, D = diag {Di, Dn) and D' = 
diag {D[, . . . , D'^) are the diagonal matrices of real eigenvalues Di and D'^. Finally, J(-D) and 
J'{D') are Jacobians resulting from the change of integration variables and the condition 
EV(Vr) < 1 translates to A < 0. 

We next 'Wick rotate' the integration over the eigenvalues from the real to the imag- 
inary axis. Since tr In (1 — ZZ) = ~ e"^^^), the Dj's can safely be sent from the 
negative real to the imaginary axis without crossing singularities. (After the rotation the 
eigenvalues of ZZ are given by e^*^\ I.e. their modulus does not exceed unity.) As a result 
of this manipulation, the integration over GL {n) /U {n) becomes one over U {n) and vice 
versa. More precisely. 



Wick 



/ f{Z, Z)= f dW f dVfiWV, V-^W) 

JMb J J 

GL(n)/U{n) U(n) 
EV(W)<1 

— > jdW J dVf{WV,V-^W) = J f{Z,Z). 



U(n) GL (n)/U (n) 
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The last line defines the new integration manifold TVb- Notice that, on the new manifold, 
the connection between Z and Z is no longer provided by Hermitian adjunction. Rather, 
we have 

Z = exp{iH) exp{H'), Z = exp{-H') exp{iH) = ZK 



APPENDIX B: ANALYTIC CONTINUATION OF THE ENERGY ARGUMENTS 

This Appendix is devoted to a discussion of the convergence problems that arise after the 
colour-flavour transformation in connection with the Gaussian integration over the bosonic 
field variables ipB- Since we will be exclusively concerned with the bosonic sector of the 
theory, all subscripts 'B' will be omitted throughout. 

That convergence problems occur follows simply from the fact that the integration over 
Z extends over field configurations with finite imaginary eigenvalues. Thus, it looks like the 
Gaussian integration over 

does not converge for all values of Zij implying that the orders of integration over ip and Z 
cannot be exchanged. 

Below we are going to show that this difficulty can be surmounted by giving the energy 
arguments Zm a finite imaginary part. More precisely, it is necessary to analytically continue 
from Zm = to = Cm ± 2z(i. Since no singularities [zm G real axis) are crossed, 

the continuation manipulation is permitted and, after it has been redone on the level of the 
subsequent Z integration, does not affect the value of the functional integral. 

To understand that a shift into the deep complex region \zm\ > 2d solves the convergence 
difficulties, we consider bilinear forms like (f)^Z(l) where Z is meant to be the BB-block of 
any matrix Zjj, 0^ = ipls^^"^ and = s^^'^ipi. To formulate the following estimates, it 
is useful to represent scalar products like the one above in terms of a bra-ket notation, 
e.g. = {(f)\Z\(f)). Now, for any normalized vector |n), {n\ZZ'^\n) < 1, a result that 

straightforwardly follows from the fact that the matrix ZZ'^ is unitarily diagonalisable and 
has (see below Eq. ([T^) eigenvalues smaller than one. Using this estimate and setting 
10) = (0|0)^/^|n), we next obtain 

\4s'i'zs'''A?^mm? = 

= {h\Z\n){n\Z^\h) (0|0) (0|0) < 

< Y.{n\Z\m){m\Z^\fL) (0|0) (0|0) = 

m 

= {h\ZZ^n) (010) (010) < (010) (010) = (^M)2, 

where {|^«)} is meant to be an extension of \n) to an orthonormalised basis. This result 
implies that, 

|Im V-Js^/^Zs^/Vil = V'J^i 
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and further 

On the other hand, after the analytic continuation of the energy arguments, 

|Im V'js^/^zs^/^Vil > 2dV»Vi- 

Consequently, the sum of the energy and the Z-contribution to the bosonic bilinear form 
has a positive imaginary part, implying that the ^/^-integration converges. 



APPENDIX C: GRADIENT EXPANSION 

This Appendix is devoted to the explicit derivation of the three contributions to the action 
(p8|). To simplify the notation we introduce the convention that all indices G A, 

whilst . . . e 5. 



1. Fluctuation action, Sq: 

Setting z = 0, ignoring massive modes, P = 1, and substituting the parameterisation 
fi into f^OD we obtain 



S,[T] = iV^strln T., +iV^strln ^ T^' . (CI) 

We next assume the lattice to be embedded in a continuum, and that the low energy field 
configurations {Tjj} can be extended to smooth functions T with continuous support. Under 
these conditions, it makes sense to Taylor expand 



T., = T, + |9.^,T, + ^C-T. + • • ■ (C2) 



where a is the lattice spacing, di^j the directional derivative along the link i ^ j, and Tj the 
function T evaluated at the coordinate of the lattice point i. Substituting this expression 
into ( |C1| ) and taking into account the fact that J^j&n^ di^jTi = 0, we obtain 



Sfi[T] = A^^strln (2dTi) + ^ str In (2c/T~^) + 



16d 



E E Tr'^^^T, + T,c^^/r-' 1 + 0{d'). (C3) 

i jeNi j i&Nj 



Up to second order in the number of derivatives (the accuracy of the present expansion 
scheme), the first two terms cancel each other: 
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str In {2dTi) + ^ str In {2dT-^) = ^ str In (T^) - ^ str In (Tj) 

i j i 3 

= E(strln(r,)-i-^strln(T, 

^ ^ ( str In (T.) - ^ 5: str In (T, + aO.^.T, + y C,T.; 

^-^y E str(-Tr^(a,,^,T,)Tr^(9,^,T,)+Tri9L,.T,) = 

= E str(9.^,T-i9._.,T. + Tr^dU.T.) = 0, (C4) 

where the last equahty is based on an 'integration by parts'. Using the fact that Z^jeA^i ^i^j = 
23"^ is twice the Laplace operator, and that in the continuum limit I^ieA ^ ^ I the re- 
maining two terms can be brought to the form 



\ i 3 

j Str {r-'d^T + Td^T-^) = I stiidT-'dT). (C5) 



16d 



2. Frequency action, Sz 



The frequency action is straightforwardly derived by expanding the action ( pOf ) to linear 
order in z. Switching directly to a continuum notation, we obtain 



Sz[T] = S[T] - S[T] 



o(z^ N{2d-lf'^ 
^=0 ' ^ Ada"^ 



str [z(T + T 



-1 



(C6) 



3. Integration over massive modes and the action, S^a- 

In order to derive the contribution S-rs^ we proceed in a standard manner and decompose 
the full action according to 

S[P, T] ^ S[X, T] =S[T] + S[X] + S,^, [T, X] , (C7) 



where we have introduced generators of the massive modes by Pij = exp(Xjj). In (|C7| ), 
S[T] is the pure Goldstone action with which we have already dealt, S[X] is purely massive, 
and Sint[X,T] is the contribution that couples massive and massless modes. As for S[X], 
the fact that the mass of the X's is of 0{N) (see below) justifies an expansion to quadratic 
order in X. With 5'int[X, T] the situation is more complicated in that we will have to 
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expect contributions C»(X"(9T, X"((9T)2)[] where n may be arbitrarily large. In order to 
justify a truncation of this series at quadratic order, n = 2, we need to employ ^ 1 
as an expansion parameter. At first sight it seems that for A^ = 0{1), we are facing a 
serious problem; Cutting the series is not a controlled approximation, whilst attempting to 
explicitly compute contributions with n ^ 1 is clearly futile. However, as discussed in the 
main body of the text, the situation is not that bad. The structure of the final result for Sra 
will be fixed by the symmetries of the action. The worst that may happen for A^ = (9(1) is a 
change in the coupling constants (from the numerical value derived below to a renormalised 
one). 

Up to an accuracy of C(X^), the action is given by 

Sm[T] = {sl!l[X,T] - ^(5«[X,T])2) , (C8) 

where Sl^^[X,T] represents the contribution to Sint[X,T] that is of nth order in X and 

(...)^ /Me-^W(...), 

is the functional average over massive modes. The starting point of our concrete evaluation 
of these expressions is the z = action 

S[X, T] = -NJ2 str In (1 + x^e^^^^) + 

+N str In j ^ e^'^T.j j + A^ ^ str In j ^ T^^^e^^^ | . (C9) 
Expanding to second order in Xij, we obtain 

E str (X.,X,,0 - ^ E str(X.,X,,), (CIO) 
^] = ^ E (A,,X,,) + ^ E Str (X,,i?,,), 

~^Estr N E XijBijX.ij + — Y BijXijBiijXiij - Y XijBiA . (Cll) 



That there are no terms of ©(X^jT) follows from the fact that the z = action depends only 
on bilinears ZZ. 
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Here we have introduced the quantities 



S,, ^ - 1 ^ i:Tr^d,^,T,. (C12) 



To prepare the perturbative integration we introduce the notation 

S^[X] = iV str (x,,Kr^,^,X,,) , (C13) 



{i',3') 



where the kernel ^ is defined through ( piO| ). The result of the integration over X will 



be essentially determined by the properties of the inverse K. For our purposes, rather than 
computing this quantity explicitly, it suffices to realize its two main features: (a) K is short- 
ranged (decays exponentially on a scale set by the level spacing). This follows from the fact 
that is massive and 'tridiagonal' on the lattice of links: More specifically, ~ A+m^, 
where A is the lattice Laplace operator and m? a mass of (9(1). (b) is symmetric. 

Functional expectation values of the contributions S[^} [X, T] can now be computed by 
means of the contraction rules 

(str (AXij)str {A'X,,j,))x = ^i^y,i'/str {AA'), 

(str {AXijA'Xi,j,))x = ^Ki.^i'j'Sti (A)str {A'). (C14) 

Instead of working out of all possible contractions of the terms appearing in ( |C11|) explic- 
itly, it is more efficient to "pre-exclude" certain contributions by virtue of their infrared- 
irrelevancy or other reasons. Specifically, the contraction of {S^^ [X, T]Y gives rise to contri- 
butions of the type ~ sir (Tj^^di^jTiT^^di'^j'TI). Since the only gradient operator compati- 
ble with the fundamental symmetries of the model is given by ~ stT{dTdT~^), contributions 
of this type can lead to no more than a renormalization of the coupling constant in S'fi[T]. 
Moreover, the structure of the kernel Kij^iij/ implies that the coupling constant does not 
change significantly, so that the contribution above can be ignored. 

As for the contraction of S^^[X,T] the 1st, 3rd, 4th and 6th contribution to that part 
of the action vanish by supersymmetry. The 2nd and 5th term give 

S^T] = {Sl^ [X, T]) = --^ ( Y: str {Tr%_^,T,) str iTr%_^,,T^K,,,,, + 



+ J2 sti{Tr^di^jTj)stiiTr%,_,^Tj)Ki,^i,j 
j,ii'eNj / 

J2 str(T-i9,^,T,)str(T-i9,^,,T,)ir.,,, 



It is a straightforward matter to check that, by the symmetry of the kernel K, only the 
contributions with j' = j and j' = oj survive, where oj is the site opposite (with respect to 
i) to j. Focusing on the two non- vanishing contributions, we obtain 
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■^J stT{T'^dT)stT{T-^dT), (C15) 



where the constant 

C = Kij^ij — Kij^ioj (^16) 
is independent of i by translational symmetry. 



4. Boundary Action, S'b 

Thus far our derivation of the action has imphcitly assumed that the system we are 
deahng with is infinitely extended. (E.g. in deriving ( PSP we have repeatedly integrated 
by parts with no account for boundary terms.) For finite size systems, however, the action 
does include boundary terms. Even for large systems, these terms must not be neglected the 
reason being that some of them are topological and hence may affect the phenomenology of 
the system on arbitrary length scales. 

To derive the boundary action, we assume that our (hypercubic) system has Ni,i = 
x,y, . . . lattice sites in x,y, . . . direction (c.f. fig. ^.) As before, the lattice spacing is set 



by a. As usual with field theories of disordered systems |^9[ the constraint that no current 
flow through the interfaces translates to the condition d±T = 0, where d± is the 

boundary 

derivative normal to the boundary. Under these conditions, the only source for potential 
boundary contributions are the first two terms in Eq. (|C^)p| As shown in Eq. (pi]), these two 
terms do not contribute to the bulk action. In order to explore their role at the boundaries, 
we imagine that the sites of the lattice are paired to doublets as indicated in fig. ^ 
Considering the case of a l(i-chain first we thus obtain, 

sI'^=^\t] = A^^strln (2rfTi) + A^^strln (2dTr^ 



1 — i—)^'^ 

N J2 (str In (T,) - str In {Tr^\)) + N ^^str In {T^,] 

i=l,3,...,Ri ^ 

.1 /-^ . „ , _ 



-N- f dxd^stT\n(T(x)) + N ^-^strln (T((Ni - l)a)) 

2 Jo 2 



N 



- (strln(T(0)) - (-)^^str In (T(L))) , (C17) 



^For reasons of operator relevancy a boundary operator must not contain more than one derivative. 
Due to current conservation, normal derivatives are not permitted. It is a straightforward matter 
to check that under these conditions the first two contributions to ( p3|) are the only ones that 
qualify as potential sources of boundary terms. 
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where i?i = A^i — 1 for A^i even and Ri = Ni — 2 for A^i odd. That scheme can straightfor- 
wardly be generahzed to higher dimensions. E.g. in = 2, we obtain 

Sf'^m = ^((-)^^+^^strln(T(Li,L2)) - (-)^^str In (r(Li, 0)) - 
-(-)^2str In (T(0, La)) + str In (T(0, 0)) 
The obvious generahsation to arbitrary dimensions is given by 



Sy>[T] = 7r, E (-)^-^^^+')^'strln(T(.iLi,...,.,L,)). 

^ Si=0,l 



(C18) 



Summing the contributions ( |C5[ ) , ( |(J6D , ( |C15| ) and ( |C18D we obtain the full continuum action 
Eq. (HD. 



APPENDIX D: THE ZERO MODE INTEGRAL 



The computation of the integral ( |32D can be most readily performed by employing the 
'polar' parameterisation of the supermatrices, 

T = kak~^, 

k = exp(^^'^^, a = exp 1^ ^ ^ , (Dl) 

where rj and /z are Grassmann variables and x,y E TZ. The invariant measure associated 
with (|D1|) is given by 

dTf{T) ^ J dr]d^ dx dy—^^^j^f{r],^,x,y). (D2) 
With this definition, the DoS takes the form 



P(e) = Po 



-oo JO smh [[x — ly) 12) 



where s+ = 7rpoe+, and 

f{x,y) 



e-"+'^ iVsitesOdd 

1 ^^sites even. 



For a lattice with an even number of sites, one can show 



1 roo /•27r 

dsP = po-Rei dx c/?/(coshxcos|/- l)e^"+(™^'^^-™^J') 

2 JO Jo 



Then, making use of the identities. 



JJs) = — / e'('s"^s/-'^?/) = _ / rfxcosh(nx) sinfscosha; - 7m/2) 
2n Jo Tc Jo 
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we find dsP = poT^^{Jo{s) — Jf(s))/2. Finally, making use of the identity 



a 



s (Jl+iis) + 4{s))] = {2n + 1) {Jl{s) - Jl_,,{s)) 



one obtains the standard result from ChRMT shown in Eq. (3^). 
For the lattice with an odd number of sites, we have the identity 



-Re 



-x+iy 



+ 



sinh {{x — iy)/2) sinh {{x + iy)/2) 



1 + 



cosh X cos y — 1 
(cosh X — cosy) 2 



As a result, the expression for the DoS separates into two contributions p(e) = Peven(e)+'5p(e)- 
The first reproduces the expression for the DoS with an even number of lattice sites 
while the second generates the correction. 



Spie) 



2n 



dx / (i|/(cosha; — cos?/)e 



js-l- (cosh x— cos y) 



Making use of the identities above we find Sp{s) = — vr^ Jo(s) Ji(s) which, taken together 
gives the expression for the DoS shown in Eq. (|35|) for the case = 1. For arbitrary A^, a 
calculation analogous to the one above obtains the more general formula (^). 
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FIGURES 




FIG. 1. The two-dimensional random flux model. 




FIG. 2. The generalized random flux model. 
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FIG. 3. Dashed: generic path contributing to @. Solid: Self retracing path which survives 
phase averaging. 



32 



FIG. 4. Visualisation of the different coupling of the C/'s and the Z's, respectively, to the field 
variables tp (represented by the solid dots). 
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FIG. 5. One- and two-dimensional finite size lattice with A'^i=odd sites in 1-direction and 
A'^2=even sites in 2-direction {2d case). The dashed ovals indicate the pairing of lattice sites 
employed in the computation of the boundary action. 
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